Large difference in the elastic properties of fee and hep hard-sphere crystals 
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We report a numerical calculation of the elastic constants of the fee and hep crystal phases of 
monodisperse hard-sphere colloids. Surprisingly, some of these elastic constants are very different 
(up to 20%), even though the free energy, pressure and bulk compressibility of the two crystal 
structures are very nearly equal. As a consequence, a moderate deformation of a hard-sphere 
crystal may make the hep phase more stable than the fee phase. This finding has implications for 
the design of patterned templates to grow colloidal hep crystals. We also find that, below close 
packing, there is a small, but significant, difference between the distances between hexagonal layers 
(c/a ratios) of fee and hep crystals. 



The simplest regular close-packed structures of hard, 
spherical particles are the face-centered cubic (fee) and 
hexagonal close-packed (hep) structures (see Fig. 0. 
Close to melting, the Hclmholtz free energies of these 
two crystal structures differ by less than 0.05% 0,0,01. 
As a consequence, hard-sphere colloids (the experimen- 
tal realization of elastic hard spheres) rarely crystallize 
directly into the more stable fee structure. Rather, crys- 
tallization initially results in the formation of a randomly 
stacked crystal [j, Q • The latter then slowly transforms 
to the stable fee structure However, pure hep 

crystals have recently been grown by colloidal epitaxy on 
patterned templates @ . At a given density, not only the 
free energies, but also the pressures and compressibilities 
of the fee and hep phases are very similar. One might 
therefore be tempted to suppose that these two crystal 
phases are similar in all their thermodynamic properties. 
Surprisingly, this is not the case. In this Letter we present 
calculations of the elastic constants of fee and hep hard- 
sphere crystals. We show that some of these elastic con- 
stants may differ by as much as 20%. As a consequence, 
a moderate deformation of the hard-sphere crystal may 
change the relative stability of the two crystal phases. 

A homogeneous deformation of a solid can be described 
by the transformation matrix a,» that relates the carte- 




o w chd 



sian coordinates Xj of a point in the undeformed solid, 
to the coordinates x\ in the deformed solid: 
where summation of repeated indices is implied. The 
(isothermal) elastic constants of a crystal are most easily 
defined in terms of an expansion of the Helmholtz free 
energy F(N, V, T) in powers of the Lagrangian strains 



FIG. 1: Sketch of the structures of the regular close-packed 
fee crystal (left) and hep crystal (right). Inequivalent close- 
packed layers are labeled with the letters A,B and, in the 
case of fee, C. An fee crystal has ABC ABC . . . stacking and 
an hep crystal has ABABAB . . . stacking. The c/a ratio is 
the distance between two close-packed layers divided by the 
distance between neighboring particles in a close-packed layer. 
The figures show the definitions of the x,y and z directions 
referred to in the text. 
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The Lagrangian strain parameters 77^ are related to the 
deformation matrix through Tjy = ^ayay — Sij). In 
Eq. the coefficients Xij(O) are simply the components 
of the stress tensor at zero deformation, Ciju are the 
second-order elastic constants, Cijkimn are the third- 
order elastic constants, and so on. For a system under 



The fee lattice has only 3 independent elastic constants 
(Cmi = Cii) C1122 = C12 and C2323 = C44 in the 
coordinate frame of the cubic unit cell). In what follows, 
we use this Voigt notation (CV,) to denote the second- 
order elastic constants. 

In order to compare the elastic constants of the fee 
and hep crystals, we used the coordinate system shown 
in Fig. ^ with the x and y directions in the hexagonal 
planes and the z direction perpendicular to these planes. 
For hep (with hexagonal symmetry) , there arc six distinct 
elastic constants, five of which are independent [lOj. To 
make a term-by-term comparison of the fee and hep elas- 
tic constant, it is convenient to ignore the full symmetry 
of the fee crystal, and only use the fact that the crystal 
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also has a lower rhombohedral symmetry. If the symme- 
try were really rhombohedral, the fee crystal would have 
six independent elastic constants. But, if we take the full 
fee symmetry into account, only three are linearly inde- 
pendent; the usual fee elastic constants can be expressed 
as linear combinations of the rhombohedral elastic con- 
stants Cy: C n = 4C' n - 3C7 33 , Ci 2 = C' 3Z + C 12 - C' n 
andC 44 = ^ 3 -i(C( 1 +Ci 2 ). 

We computed the elastic constants by calculating the 
stress response to a small applied strain, using molecu- 
lar dynamics simulations [llj. At zero deformation, the 
stress response of a system with isotropic pressure P is 
given by a generalization of Hooke's law: 
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{SijSki - SuSjk - SjiS ik )P + Cijki (2) 



For the MD simulations, we used the event-based algo- 
rithm described by Rapaport0l- The pressure tensor is 
calculated as the time average of the dyadic product of 
the collisional momentum exchange vector and the par- 
ticle separation vector for each two-particle collision [l3|. 

We performed simulations for a range of amplitudes of 
each type of deformation. The second-order elastic con- 
stants were deduced from the linear part of the stress- 
strain relation. In principle, all elastic constants can 
also be calculated in a single simulation using fluctua- 
tion methods 0, ITEl Put. However, these methods suffer 
from slow convergence [l4|. We found the stress-strain 
method to be the most efficient. 

For some deformations, we also computed the third- 
order elastic constants from the second derivative of the 
stress tensor with respect to deformation: 
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28tu$rsTij + (SitSj r + 5i r 8jt)T su 

—S sr [SjtTi u 4- 5aT U j) 

~S sr Cijtu H~ $irC sjtu ~T~ $jrCist u 
~\~S r tCijsu H~ Cij rs t u 



(3) 



The third-order elastic constants Cij rs t u appear in the 
last term. 

The simulations were performed on systems with 6 x 6 x 
6 = 216, 12x12x12 = 1728 and 24x24x24 = 13824 par- 
ticles. The maximum applied deformation at lower den- 
sities was 4 ■ 10 -3 ; higher densities required even smaller 
deformations to keep the stress response linear. The mea- 
sured elastic constants between the melting point (pack- 
ing fraction <f) — 0.54329 01) an d close packing are given 
m table H| 

At all densities, the values of the fee and hep elas- 
tic constants differ significantly (see Fig. The rela- 
tive differences between the elastic constants appear to 
remain approximately constant over the entire density 
range. The largest difference between fee and hep (up 
to 20%) was found for C 12 - Yet, the compressibilities of 
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FIG. 2: Relative difference A = \C'^ C - C% cp \/C'£ c between 
fee and hep C'n, C[ 2 an d C' 13 elastic constants as a function 
of packing fraction cj>. The results shown were obtained in 
simulations of a system of 216 particles, with c/a = ^/8/3. 
The curves only serve as guides to the eye. 
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FIG. 3: Variation of the fee and hep free energy with de- 
formation £ (see Eq. |IJ . The symbols indicate the results of 
Einstein free energy calculations (see text). The continuous 
curves were computed on basis of the calculated second and 
third order elastic constants. The error bars of the Einstein 
free-energy calculations are about one sixth the size of the 
symbols. The two horizontal arrows show the predictions for 
A/ obtained by neglecting the third-order elastic constants. 



the two phases are identical to within the measurement 
error. For instance, at melting: Kj cc = 0.02422(5) vs. 
K% cp = 0.02424(5) (for 1728 particles). We computed 
these compressibilities in two ways: (a) from the appro- 
priate linear combination of elastic constants and (b) di- 
rectly from the equation of state 0] . The results are the 
same, to within the statistical error. At the same den- 
sity, the pressures of the fee and hep phases are also very 
similar: P fcc = 11.568(1) and P hcp = 11.571(1). Finally, 
the free energies differ only by about 1.12(4) lO~ 3 ksT 



per particle [lj, 0, ■ 

The difference between the fee and hep elastic con- 
stants is surprising because, already in 1967, Stillinger 
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TABLE I: Second-order elastic constants of fee and hep hard-sphere crystals at densities between the melting point and close 
packing. The values for the hep structure with c/a = are shown in upright font. The (almost identical) results for a 

fully relaxed c/a ratio: c/a = y/8/3(l - 7.5 • 10~ 4 ) at <f> = 0.543, are shown in italics. The simulation equilibration time was 
1 • 10 4 collisions per particle. Data were collected during typically 2 • 10 6 collisions per particle for the 216 particle system, and 
6 ■ 10 4 collisions per particle for the 13292 particle system. For each deformation, 8 simulations were done at different strain 
amplitudes to check linearity of the stress response. The calculations of the stress-strain curve for each type of deformation 
involved simulations totaling several billion collisions 6.4 • 10 9 collisions (one week on an Athlon 1600+ CPU). 
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TABLE II: Values for the computed third-order elastic con- 
stants at melting (<j> — 0.54329). These numbers were ob- 
tained for fee and hep systems containing 13292 particles. 
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and Salsburg [l9( had pointed out that a simple free- 
volume model predicts that the fee and hep elastic con- 
stants should be equal. However, they also showed that 
pair and triplet correlation effects can lead to differences. 
Still, we were surprised by the magnitude of the com- 
puted differences, in particular for C[ 2 . To double-check 
our calculations of the elastic constants, we performed a 
second, fully independent calculation where we directly 
computed the free energy of the crystals in various states 
of deformation. The free energy of the (deformed and un- 
deformed) crystals was calculated using a 20-point Ein- 
stein integration 01 ■ We found that the results obtained 
by the two methods were completely consistent. For ex- 
ample, in Fig. 03 we show the results of the two calcula- 
tions for free energy change due to a deformation of the 
form 

/l+C 0\ 

1/(1 + (4) 
V 1/ 

To lowest order in f, AF/V = (~2T XX + Cn-C 12 )£ 2 , for 
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FIG. 4: Equilibrium anisotropy {v = 1 — £ /y§) for the hep 
hard sphere crystal as a function of packing fraction. 



this deformation. As the figure shows, the differences in 
elastic constants C' xl and C[ 2 , for fee and hep, are so large 
that a deformation of 1.2% is enough to make hep more 
stable than fee. The free energy increase of the fee phase 
due to a deformation of 2% is Af fcc = 1.93(1) ■ 1(T 2 , 
while for hep it is only Af hcp = 1.66(1) • 1(T 2 . Fig. 
also shows the effect of the third-order elastic constants. 
To within the statistical accuracy of our simulations, the 
relevant third-order elastic constants (see table ITT|> . were 
found to be the same for fee and hep. Hence, they do not 
affect the free energy difference between the two lattices. 

For the undeformed fee system, all three diagonal com- 
ponents of the pressure tensor are equal. However, this 
does not hold for a hep system at the same c/a-ratio 
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(i.e. for the same spacing between the close packed [111]- 
planes). If we fix the c/a ratio at the fee value (-^/8/3), 
the stresses exhibit a slight anisotropy. For the 13292 
particle system T xx and T yy are equal (as they should): 
T xx = -11.587(1), Tyy = -11.588(1). However, T zz is 
significantly different: T zz — —11.537(1). FromEq.^ we 
can derive what change in the c/a ratio is needed to make 
the pressure isotropic. We find that, at melting, isotropy 
is restored for a c/a ratio of y/8/S(l - 7.5(2) • 10" 4 ). At 
higher densities, this value approaches the close-packing 
value c/a = \/8/3, as can be seen in Fig.Q] Stillinger and 
Salsburg used the cell-cluster method to estimate the 
difference of the fee and hep c/a ratios. Our simulations 
show that, close to melting, the effect is one order of 
magnitude larger than predicted. The free energy differ- 
ence between the equilibrium hep and fee crystals is only 
slightly changed by this relaxation of the hep c/a ratio: 
it becomes 1.050(5) • 10" 3 fc s T per particle for N = 13292 
at melting. 

As can be seen from the results in table |U for <j) = 
0.543 — where the c/a ratio differs most from fee — the 
effect of relaxing c/a to its equilibrium value, is barely 
significant. For this reason, most hep elastic constants 
in table |U were computed for c/a = y/8/3. The table 



also shows that the elastic constants depend somewhat 
on system size, but the effect is too small to change the 
qualitative picture. 

In colloidal-epitaxy experiments [9j , the best hep crys- 
tals were obtained when the patterned template was 
stretched by 2.6% with respect to the expected lattice 
spacing at the experimental packing fraction ((f> = 0.68). 
The templates used matched a diagonal cut through the 
xy plane of Fig. ^ Together with the stress produced by 
gravity (resulting in a strain perpendicular to the tem- 
plate plane of —2.8%), this strain is comparable to the 
strain of Eq.^and would result in a free energy difference 
of about 3 ■ 10 _2 fc bT per particle in favor of hep. The 
present simulation results may help experimentalists in 
designing optimal templates to grow selectively colloidal 
hep or fee crystals. 
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